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ABSTRACT 

Gamma-rays from the decay of A1 offer a stringent constraint on the Galaxy’s global 
star formation rate over the past million years, snpplementing other methods for qnanti- 
fying the recent Galactic star formation rate, such as equivalent widths of Ha emission. 
Advantages and disadvantages of using ^®A1 gamma-ray measurements as a tracer of the 
massive star formation rate are analyzed. Estimates of the Galactic ^®A1 mass derived 
from COMPTEL measurements are coupled with a simple, analytical model of the ^®A1 
injection rate from massive stars and restrict the Galaxy’s recent star formation rate to 
5 ± 4 Mq yr“h In addition, we show that the derived ^®A1 mass implies a present day 
Type II -|- Ib supernovae rate of 3.4 ± 2.8 per century, which seems consistent with other 
independent estimates of the Galactic core collapse supernova rate. If some independent 
measure of the massive star initial mass function or star formation rate or Type II -|- Ib 
supernovae rate were to become available (perhaps through estimates of the Galactic ®°Ee 
mass), then a convenient way to restrain, or possibly determine, the other parameters is 
presented. 


Subject headings: gamma rays: theory - nucleosynthesis, abundances - 
supernovae: general 
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1. INTRODUCTION 

Through its 1.809 MeV gamma-ray line, Galactic was discovered in 1979 with 
the High Energy Astronomy Observatory-C spectrometer (Mahoney et ah 1982). Sev¬ 
eral measurements of the integrated 1.809 MeV flux have been performed since then, as 
reviewed by Prantzos &: Diehl (1996). The most reliable of these measurements are prob¬ 
ably derived from data obtained from the Gamma-Ray Spectrometer aboard the Solar 
Maximum Mission spacecraft (Share et al. 1985; Harris, Share, &: Leising 1994), and the 
COMPTEL Imaging Telescope aboard the Compton Gamma-Ray Observatory (Diehl et 
al. 1995a). All estimates of the absolute ^®A1 mass in the Galaxy rest on assumptions 
about the spatial distribution of the sources, as the 1.809 MeV measurements themselves 
do not carry distance information. This situation may change with future high-resolution, 
high-sensitivity instruments if the line shape and Doppler shift of the 1.809 MeV line can 
be extracted (Gehrels & Ghen 1996). 

Erom the HEAO-G data, Mahoney et al. determined 3 Mq of ^®A1, assuming the smooth 
spatial distribution derived from COS-B measurements of Galactic gamma-rays in the 100 
MeV regime. All other non-imaging instruments basically confirm this result, based on the 
same (or equivalent) assumptions about spatial source distributions. It became evident, 
however, with the COMPTEL imaging data (~ 5° degree spatial resolution) that the dis¬ 
tribution of 1.809 MeV emission is significantly different than that derived from COS-B 
measurements (Diehl et al. 1995a,b). Although the ridge of the Galactic plane dominates 
the emission, several prominent regions of emission (particularly Cygnus and Vela), sug¬ 
gest substantial deviations from a smooth emission pattern following the distribution of 
gas. The GOMPTEL team fitted several models of candidate source distributions to the 
COMPTEL data, such as CO survey data (Dame et al. 1987), analytical models based on 
exponential disks, and H H region data evaluated in the context of spiral-arm structure 
(Taylor & Cordes 1993). If the prominent, localized regions of emission beyond the inner 
Galaxy are excluded from a fit to such models, then all fits to axisymmetric models yield 
a Galactic mass of 3 ± 0.5 Mq (Diehl et al. 1995; Knodlseder et al. 1996). However, 
the marked asymmetry in the emission profile along the disk suggests that spiral struc¬ 
ture is important (Prantzos & Diehl 1996; Chen, Gehrels, &: Diehl 1995). If one adopts a 
composite model of disk-like emission plus contributions from sources along spiral-arms, 
then fits to the GOMPTEL data results in a total Galactic mass of 2.7 ± 0.8 Mq, with 
0.7 ± 0.3 Mq attributed to the spiral arm component (Knodlseder et al. 1996). 

Candidates for the origin of ^®A1 include massive stars (through their supernovae ejecta 
and Wolf-Rayet wind phase contributions), asymptotic giant branch (AGB) stars, and 
classical novae with metal enriched atmospheres (see Prantzos & Diehl 1996 for details) 
and While yields from all these candidate Al sources are uncertain, it seems safe to assume 
that core-collapse supernovae as the dominant source will not be challenged. Thus, the 
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COMPTEL measurements and analysis can be interpreted as: 0.7 Mg of ^®A1 if only the 
spiral-arm component is assigned to Type II -|- Ib supernovae, to 2.5 Mg of ^®A1 if ~ 80% 
of the emission is assigned to Type Il/Ib supernovae and no foreground contributions from 
localized emission regions similar to Vela or Cygnus lie in the direction of the inner Galaxy. 
Since the decay time of ^®A1 (ri /2 = 7.5 x 10® yr) is short compared to Galactic rotation 
timescales (rgai ~ 10 ® yr) this estimated ^®A 1 mass range serves as an important constraint 
of the stellar population responsible (i.e., massive stars) for the synthesis. 

One of the recurring concepts of this paper is the relationship between the derived A1 
mass estimates, initial mass function (henceforth IMF), recent global star formation rate 
(henceforth SFR), and present epoch Type II -|- Ib supernovae rate. The connections made 
in this paper between these four quantities are new, or at least not widely recognized, but 
the raw science of the four quantities themselves is not new. It is worthwhile in making the 
connections to succinctly summarize the shape of the IMF, tracers of the Galactic SFR, 
and stellar core-collapse estimates. 

2. STAR FORMATION RATE GONSTRAINTS 
2.1. The Initial Mass Function 

The hrst empirical determination of the observed IMF showed that the number of stars 
between 0.4 Mg and 10 Mg could be described as a power-law with the index 7 = —2.35 
(Salpeter 1955). Studies since then (Miller & Scalo 1979; Humphreys & McElroy 1984; 
Scalo 1986; Rana 1991; Tinney, Mould & Reid 1992; Parker & Garmany 1993; Reid 1994; 
Hunter 1995; Massey et al. 1995a,b; Kroupa 1995; Mera, Ghabrier & Baraffe 1996; Mayya 
& Prabhu 1996; Hunter et al. 1996) suggest that the observed IMF becomes flatter than 
a pure power-law at the smallest stellar masses (7 ^ — 1 for m < 0.5 Mg) and becomes 
steeper for the most massive stars (7 ^ —3.3 for m > 10 Mg). Some studies indicate that 
the IMF has more structure than either a power-law or log-normal form (Rana 1991), 
while others argue that the IMF is closer to a power-law and hence has less structure 
(Scalo 1986). Overall, the shape of the IMF appears to be quite robust (centered on the 
Salpeter —2.35 exponent) and seems not to change very much from one star-forming region 
to another. In any event, proceeding from the observed luminosity function to the implied 
IMF depends upon the stellar evolutionary tracks used in the htting procedure (Tinsley 
1980; Elmgreen 1995a,b; Efremov 1995; Adams & Fatuzzo 1996; Arnett 1996). 

2 . 2 . Galactic Star Formation Rates 

Despite many uncertainties, the results of various studies (Schmidt 1959, 1963; Searle, 
Sargent & Bagnuolo 1973; Larson & Tinsley 1974; Gohen 1976; Huchra 1977; Smith, 
Biermann &: Mezger 1978; Lequeux 1979; Talbot 1980; Tinsley 1980; Kennicut 1983; 
Kennicut & Kent 1983; Glisten & Mezger 1983; Turner 1984; Lacey & Fall 1985; Dopita 
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1985; Gallagher, Bushouse, & Hunter 1989; Romanishin 1990; Rana 1991; Kennicut 1992; 
Lada 1992; Gallagher & Gibson 1993; Kennicut, Tamblyn &: Congdon 1994; Hill, Madore, 
& Freedman 1994; Gallager & Scowen 1995; Gallego et ah 1995; Lada & Lada 1995) lead 
to the picture of ^ 10% of the current SFR occuring in the innermost 1 kpc of the Galaxy, 
and most of the remaining 90% concentrated between 5 and 9 kpc from the center, which 
is where most of the Galaxy’s giant molecular clouds, infrared emission, and other signs of 
intense star formation reside. The current SFR for the whole Galaxy has been estimated 
to be 0.8 Mq yr^^ (Talbot 1980), 3.0 Mq yr^^ (Turner 1984), 5.3 Mg yr^^ (Smith et al. 
1978), 13.0 Mg yr“^ (Giisten h Mezger 1982), and 6.0 Mg yr“^ (Pagel 1994). Glisten h 
Mezger (1982) also estimated the massive star SFR in the spiral-arms to be 5 ± 2 Mg yr^L 
Note that all of these estimates use various indicators of the massive star population, and 
then convert these indicators to a total SFR by means of an assumed (universal) IMF. 

2.3. Ha Line Widths 

Equivalent widths of Ha emission have been the best available and most popular 
method for quantifying the present SFR. This is because Ha equivalent widths are di¬ 
rectly proportional to the number of Lyman continuum photons emitted by massive stars 
and hence proportional to the SFR. Other SFR measures or indicators (H/5, Hy, [O HI] 
A5007, [O H] A3727, integrated UBV colors, infrared luminosities, IRAS fluxes, free-free 
radio emission from H H regions, magnetic held strengths, and brightest individual star 
counts) are more affected by stellar absorption, interstellar reddening, excitation strength, 
metallicity, dust abundances, dust composition, incompleteness, sky coverage, or resolu¬ 
tion limitations than Ha emission (see references above). Even when Ha measurements 
are combined with some of the alternative indicators, the derived SER is a lower limit since 
even Ha is not completely immune (just less sensitive) from the contaminants listed above. 
It has been suggested that near-infrared recombination lines of Bry could be an even better 
measure of the current SER (Leitherer & Heckmann 1995), but instrumentation difficulties 
impede progress along this avenue at present. 

2.4. Gamma-ray Measurements 

Gamma-rays from the decay of ^®A1 offer a unique measure of the present SER in the 
Galaxy. Several of the difficulties noted above are mitigated by the transparency of the 
Galaxy to gamma-rays (e.g., absence of interstellar reddening), but several difficulties re¬ 
main (e.g., spatial resolution limitations). Nevertheless, gamma-rays offer a complimentary 
indicator of the Galaxy’s present epoch SER. 
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The IMF by number (assumed universal and constant independent), the normalization 
condition, and the normalization constant are 
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( 1 ) 


respectively. The mean mass of stars, the fraction of all stars which become core-collapse 
supernovae (Type II -|- Type Ib), and the steady-state core-collapse supernovae rate, re¬ 
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where is the stellar birthrate in number per year, T is the SFR in Mq per year. Ml is the 
smallest stellar mass in the distribution. Mu is the largest stellar mass in the distribution, 
and Msn is the smallest stellar mass which undergoes core-collapse. The mean yield of 
^®A1 from Type II -|- Ib events is 
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Finally, the steady state injection rate of ^®A1 is 


M26 — R-sn < y > — 


^ Yeff 

< m > 


[ Me 1 

1 year J 


( 3 ) 


( 4 ) 


Eq. (4) may be solved for the global SFR, T, for a given observed ^®A1 mass in 
a steady state galaxy, and a given IMF exponent 7 . The results of such a procedure 
is shown in the lower panel of Figure 1. Each labeled curve corresponds to a different 
Galactic ^®A1 mass (in solar masses), with the preferred ^®A1 mass range (0.7 - 2.8 Mg) 
imposed by the COMPTEL observations shown as the grey band. Integration limits of 
Ml = 0.1 Mg, Mu = 40 Mg, and MgN = 10 Mg were used in constructing Eig. 1, but the 
chief conclusions are quite robust with respect to reasonable variations in the integration 
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limits. The mean and effective ^®A1 yields in eqs. (3) and (4) were calcnlated with the 
^®A1 mass ejected in the Woosley & Weaver (1995) massive star models. There is about an 
order of magnitude difference in the A1 yields if the results of the Thielemann, Nomoto, & 
Hashimoto (1996) survey are used instead of Woosley &: Weaver. The bulk of the synthesis 
of this radioactive isotope takes place in the presupernova star. It is imperative to follow 
this stage of the star’s evolution with a sufficient nuclear reaction network, especially 
during the last few hours of convective neon and oxygen burning. Woosley & Weaver used 
a 200 isotope network from the main sequence through the explosion, while Thielemann 
et al. follow the presupernova evolution from an initial helium core mass with an a-chain 
network. Only during the explosive phases of the evolution do Thielemann et al. switch to 
a larger reaction network. This accounts for most of the difference in the ^®A1 production 
in the two surveys. Deviation from straight lines in the lower panel of Fig. 1 is due to 
the IMF exponent approaching the removable singularity at 7 = — 1 (see eq. 1). Only a 
mathematical reason, not a physical one, is responsible for the flattening of the curves. 

The horizontal dimension of the dashed box in the lower panel of Fig. 1 is centered 
on the Salpeter —2.35 exponent, and is representative of the range of IMF exponents for 
massive stars encountered in the literature. Vertical dimensions of the dashed box were set 
by requiring consistency between the COMPTEL estimates of the Galactic ^®A1 mass and 
the simple model for the Al injection rate from massive stars. The lower panel of Fig. 1 
suggests that the global SFR in the Galaxy during the past million years is restricted to 
5 ± 4 Mq yr“h This is consistent with the Giisten & Mezger (1982) Ha estimate of the 
massive star SFR in the spiral arms of 5 ± 2 Mg yr^^, and the more recent determinations 
of the Galaxy’s global SFR (see §1). 

Eq. (2) may be solved for the core-collapse supernova rate given the global SER 
4/ and the IMF exponent 7 . This solution is shown in the upper panel of Fig. 1, for 
the SFRs calculated in constructing the lower panel of Fig. 1. As before each labeled 
curve corresponds to a different Galactic ^®A1 mass, with the preferred ^®A1 mass range 
(0.7 — 2.8 Mg) imposed by the GOMPTEL observations shown as the grey band. Here 
the curves are straight lines (as expected) since the approach to the removable singularity 
at 7 = —1 is embedded in both factors (< m >, Esn) of eq. ( 2 ) and they cancel each 
other. Eor the same plausible range of IME exponents considered above, the GOMPTEL 
estimates of the Galactic ^®A1 mass appear to imply a core-collapse supernovae rate of 
3.4 ± 2.8 per century. 


3. DISGUSSION 

Direct measurement of the Galactic supernova rate is difficult owing to possible incom¬ 
pleteness in historical observations, and uncertainty as to the fraction of the Galactic disk 
and altitude that are sampled. Indirect inference from supernova rates in similar galaxies 
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is adversely affected by the imprecise value of the Hubble constant, and the uncertainty in 
estimating the total blue luminosity and morphological classihcation of our Galaxy. Sys¬ 
tematic searches for extragalactic supernova are also hampered by the need to know the 
distance, luminosity, and Hubble class of the host galaxy, as well as the dates and limiting 
magnitude of each observation. Such detailed information is only available in a few dozen 
supernova catalogs. Based on these surveys, estimates of the core-collapse and thermonu¬ 
clear driven supernova rates were derived and discussed by van den Bergh &: Tammann 
(1991), and Cappellaro (1993). These estimates assumed that the peak luminosity of each 
supernova class was a standard candle, and a large correction for edge-on spirals (sini ef¬ 
fect). Using the extragalactic estimates with a total Galactic blue luminosity of 2.3 x 10^° 
Lq, a Hubble constant of 75 km s“^ Mpc“\ and a Sbc Galactic morphology, the Galactic 
core-collapse supernova rate has been estimated to be 4.1 per century (van den Bergh &: 
Tammann 1991) and 2.4 — 2.7 per century (van den Bergh & McClure 1994; Tamman, 
Lofiler, &: Schroder 1994). These estimates agree (perhaps auspiciously) with the core¬ 
collapse supernovae rate implied by a near-Salpeter IMF exponent and the COMPTEL 
derived ^®A1 mass. 

While the general agreement found between estimates of the COMPTEL derived ^®A1 
mass, the range of massive star IME exponents encountered in the literature, complimen¬ 
tary measures of the recent SFR, and the present epoch Type H -|- Ib supernovae rate 
may be fortuitous and reminiscent of epicycles, it does point to a consistent picture. Less 
speculative is the fact that gamma-rays from the decay of certain radioactive nuclei, such 
as ^®A1 and ®°Fe, offer a unique measure of the present SFR in the Galaxy that is compli¬ 
mentary to other popular indicators of the Galaxy’s present epoch SFR (e.g, Ha, H/3, Hy, 
[O HI] A5007, [O H] A3727, integrated UBV colors, infrared and radio luminosities, and 
stellar counts). If some independent measure of the massive star IMF exponent or SFR or 
Type H -|- Ib supernovae rate were to become available (perhaps through measurements 
of the Galactic ®°Fe mass), then Fig. 1 offers a convenient way to constrain, or possibly 
even determine, the other parameters. 
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FIGURE AND CAPTION 
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Fig. 1.— Global star formation rate (lower panel) and current Galactic Type II + Ib 
supernova rate (upper panel) vs the IMF exponent. Each labeled curve corresponds 
to a different Galactic ^®A1 mass (in solar masses), with the preferred ^®A1 mass range 
(0.7—2.8 Mq) imposed by the COMPTEL observations shown as the grey bands. The 
horizontal dimension of the dashed boxes are centered on a Salpeter —2.35 exponent, 
and are representative of the range of IMF exponents for massive stars encountered 
in the literature. Vertical dimensions of the dashed boxes were set by requiring con¬ 
sistency between the COMPTEL estimates of the Galactic ^®A1 mass and the simple 
model for the ^®A1 injection rate from massive stars. This consistency then appears 
to imply a Galactic SFR during the past million years of 5 ± 4 Mq yr“\ and a core¬ 
collapse supernovae rate of 3.4 ± 2.8 per century. 
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